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Abstract: We propose a new family of matrix models whose 1/N expansion captures the all¬ 
genus topological string on toric Calabi-Yau threefolds. These matrix models are constructed 
from the trace class operators appearing in the quantization of the corresponding mirror curves. 
The fact that they provide a non-perturbative realization of the (standard) topological string 
follows from a recent conjecture connecting the spectral properties of these operators, to the 
enumerative invariants of the underlying Calabi-Yau threefolds. We study in detail the resulting 
matrix models for some simple geometries, like local and local F 2 , and we verify that their 
weak’t Hooft coupling expansion reproduces the topological string free energies near the conifold 
singularity. These matrix models are formally similar to those appearing in the Fermi-gas for¬ 
mulation of Chern-Simons-matter theories, and their 1/N expansion receives non-perturbative 
corrections determined by the Nekrasov-Shatashvili limit of the refined topological string. 
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1. Introduction 

One of the most surprising aspects of string theory is that, in some circumstances, it can be 
described by very simple quantum systems. For example, non-critical (super)strings can be 
formulated in terms of double-scaled matrix models or matrix quantum mechanics. In some cases, 
these equivalent descriptions provide as well a non-perturbative definition of the corresponding 
string theory models. More recently, some simple quantities in fully-fledged superstring theory, 
like partition functions, have been also expressed in terms of matrix integrals, by combining the 
AdS/CFT correspondence with supersymmetric localization. 

It has been suspected for a long time that topological strings on Calabi-Yau (CY) mani¬ 
folds should be also described by simple quantum models. It was found in [1] that the type B 
topological string on a special class of non-compact B-model geometries has such a description, 
in terms of conventional matrix models. The CY backgrounds considered in [1] are useful from 
the point of view of engineering supersymmetric gauge theories, but they have no mirror geome¬ 
tries and no enumerative content. It was later found that bona fide topological strings on An 
fibrations over can be described by Chern-Simons matrix models [2], as a consequence of the 
Gopakumar-Vafa large N duality [3] and its generalizations [4]. 

Recently, a correspondence has been proposed between topological strings on toric CY three¬ 
folds, and the spectral theory of operators arising in the quantization of their mirror curves [5]. 
The idea that the topological string free energies could emerge from the quantization of mirror 
curves was first proposed in [6]. In [7, 8], building upon the work of [9], it was shown that 
a perturbative treatment of the quantum mirror curve leads to the Nekrasov-Shatashvili (NS) 
limit of the refined topological string. However, building on the study of matrix models for 
Chern-Simons-matter theories [10] and their AdS/CFT duals [11] (see [12] for a review), it was 
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Figure 1: Given a toric Calabi-Yau threefold, the quantization of its mirror curve leads to a trace class 
operator p. The standard topological string free energy is obtained as the’t Hooft limit of its fermionic 
traces Z{N,h). 


pointed out in [13] that the standard topological string also emerges from the quantum curve, 
once non-perturbative corrections are taken into account. 

The proposal of [5] incorporates all these ingredients in an exact treatment of the quantum 
curve. According to [5], to each mirror curve of a toric Calabi-Yau threefold, one can associate 
a positive, trace class operator on L^(M). This was rigorously proved for a large number of 
geometries in [14]. Therefore, these operators have a positive, discrete spectrum, and their 
Fredholm or spectral determinants are well dehned. In [5], an explicit formula for these spectral 
determinants was conjectured, involving both the NS limit of the rehned topological string and 
the conventional topological string. The conjectural, exact formula of [5] leads in addition to 
an exact quantization condition determining the spectrum, which generalizes previous studies 
of the spectral problem [13, 15, 16, 17] and is conceptually similar to other exact quantization 
conditions appearing in Quantum Mechanics (see for example [18]). This establishes a novel and 
precise connection between spectral theory and mirror symmetry. 

One of the consequences of the correspondence of [5] is that the conventional topological 
string free energy (at all genera) appears as a ’t Hooft limit of the spectral determinant. A 
very useful way to encode the information in the spectral determinant is in terms of the so- 
called fermionic spectral traces Z{N,h) (see section 2.1 for precise dehnitions). As we will show 
in detail, these traces have a natural matrix model representation. It then follows from the 
conjecture of [5] that the’t Hooft limit of this matrix model. 


N —)• oo, h —)■ oo, 

is given by the asymptotic expansion 


— = A fixed, 
h 


( 1 . 1 ) 


iogZ(iv,h) = J]j-,(A)h2-2^ 
<?>0 


( 1 . 2 ) 


where J^g{X) are the standard topological string free energies, and the ’t Hooft parameter A is 
a flat coordinate for the CY moduli space^. Therefore, the conjecture of [5], together with the 

^In this paper, as in [5], we will focus on local del Pezzo Calabi-Yau’s, where there is a single “true” modulus, 
and correspondingly a single’t Hooft parameter. 
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representation of Z{N,h) in terms of a matrix integral, gives a matrix model for topological 
strings on toric Calabi-Yau threefolds. The construction is summarized in Fig. 1. We should 
add that the resulting matrix model is a convergent one, i.e. the matrix integral is well-defined, 
as a consequence of the operator being of trace class. 

As explained in [5], the correspondence between spectral theory and mirror symmetry pro¬ 
vides a rigorous, non-perturbative completion of the topological string, in the sense that the 
genus expansion of its free energy, which is known to be a divergent series, is realized as the 
asymptotic expansion of a well-defined function. The implementation of the correspondence of 
[5] that we are presenting here, in terms of matrix models, makes this point particularly clear; 
the quantity Z(N, K) is manifestly well-defined for any positive integer N and any real h, since it 
is defined by the spectral theory of a trace class operator (it can be also analytically continued 
to complex N, as first explained in [19], and it is likely that it has an analytic continuation to 
complex values of h). The’t Hooft expansion of Z{N, h) is given exactly by the genus expansion 
of the topological string free energy. However, there are non-perturbative corrections to the ’t 
Hooft expansion, due to large N instantons, which are also predicted by the conjecture in [5], 
and they are encoded in the NS limit of the refined topological string. 

The fact that the matrix model representing Z(N, H) leads to the topological string free 
energies is not at all obvious. On the contrary, it is a highly non-trivial prediction of the conjecture 
of [5]. The’t Hooft limit of Z{N, h) probes the strong coupling limit oi the spectral problem (since 
h is large), and in particular the non-perturbative instanton corrections to the perturbative WKB 
expansion. For this reason, in this paper we will perform detailed calculations in some examples to 
verify that, indeed, the’t Hooft expansion of the matrix model defined by the trace class operator 
gives the topological string free energies. This constitutes an analytic test of the instanton 
corrections postulated in [5]. 

As we mentioned above, the large N duality of Gopakumar-Vafa [3] and its generalizations 
[4] provide a matrix model representation of the free energies of topological strings in certain 
geometries, as it was tested in [4, 20, 21]. Although in this paper we focus on local del Pezzo 
geometries, our matrix models are potentially valid for any toric geometry, in contrast to the 
duality of [3, 4]. For example, we will study in detail a matrix model for local which has no 
counterpart in the framework of [3, 4]. It would be interesting to understand the relationship 
between the matrix models for topological strings obtained in [4] and the ones described here. 

Matrix models describing topological strings on more general backgrounds were also proposed 
in for example [22, 23, 24, 25, 26]. We should note that these models are very different from 
the ones we construct here: first of all, they are engineered ah initio to reproduce formally the 
topological string free energies; in contrast, our matrix models are defined by the trace class 
operators obtained in the quantization of the mirror curve, and the fact that they lead to the 
correct topological string free energies is a consequence of the non-trivial conjecture of [5]. Second, 
in the models of [22, 23, 24, 25, 26], the rank of the matrix N plays an auxiliary role, while in our 
case it is a flat coordinate for the Calabi-Yau, as in other large N dualities. Third, the matrix 
models in [22, 23, 24, 25, 26] are often formal (i.e. not convergent) and therefore can not define 
a non-perturbative completion of the theory; our matrix models are convergent and lead to a 
non-perturbative completion. 

This paper is organized as follows. In section 2 we review elementary aspects of trace class 
operators on L^(M) and we note that their fermionic traces have matrix model-like representa¬ 
tions. We then focus on the operators coming from quantized mirror curves, and write down 
explicit matrix models for some of them, including the ones relevant for local and local F 2 . 
These models can be studied in the ’t Hooft expansion, and we compute their weakly coupled 
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expansion to the very first orders. In section 3 we review the conjecture of [5] and we spell out 
in detail its prediction for the’t Hooft expansion of the fermionic traces. Then, we test the con¬ 
jecture in detail for the matrix models describing local and local F 2 . We conclude in section 
4 and we list some open problems for the future. In the Appendix, we list some results for the 
weakly coupled expansion of the matrix model free energies. 


2. From operators to matrix models 


2.1 General aspects 

Let us then begin with some general aspects of operator theory. Let p be a positive-definite, 
trace class operator on the Hilbert space % = L^(M), depending on a real parameter h. Since its 
eigenvalues are discrete and positive, we will denote them by e~^", where n = 0,1, • • •. Due to 
the trace class property, all the spectral traces of p 

Z, = Tr^/ = £ = 1 , 2 ,---, (2.1) 

n>0 


exist. One can also define the fermionic spectral traces as 

Z{N,h) = Tx^Mp^^{K^{p)), iV = l,2,. 


( 2 . 2 ) 


where the operator {p) is defined by p^^ acting on ifH) (see for example [27] for a review 
of these constructions)^. The fermionic spectral traces are related to the standard spectral traces 
(2.1) by the equation 

{me} e 

where the ' means that the sum is over the integers mi satisfying the constraint 

Y,£mi = N. (2.4) 


The Fredholm or spectral determinant of p can be defined as 

00 

H(k, h) = det(l -|- np) = (l + , 

n=0 

or, equivalently, by the expansion around k = 0, 

00 

E{K,h) = 1+ Z{N,h)K^, 


(2.5) 


( 2 . 6 ) 


7V=1 


and is an entire function of k [27]. Note that, if p is interpreted as a one-particle thermal density 
operator, the fermionic trace Z{N,h) is the canonical partition function for an ideal Fermi gas 
of N particles. It then follows that Z(N,h) has the matrix-model-like representation 

^ ^ J X det {p{xi, Xj)) . (2.7) 

CTSSiV * 

^The terms “spectral trace” and “fermionic spectral trace” do not seem to be standard, but since these objects 
are of paramount importance in our construction, we had to find a name for them. 
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In this equation, p{xi,X 2 ) is the kernel of the operator p, 

p(xi,X2) = (xi|p|x2), (2.8) 

Sn is the permutation group of N elements, and e{a) is the signature of a permutation a G Sn- 
The above formula encapsulates the fermionic nature of Z(N,h). 

At this stage, calling (2.7) a matrix model might seem excessive. However, the integrand 
of (2.7) has an essential property, typical of the integrand of a matrix model in the eigenvalue 
representation: it vanishes whenever Xi = Xj. In standard matrix models, this is an indication 
of eigenvalue repulsion. To further understand this parallelism, note that, when solved in terms 
of orthogonal polynomials, the partition function of a Hermitian matrix model is given by the 
expression 

^ y d^xdet{KNixi,Xj)), (2.9) 

where the kernel K]sf(x, y) can be written in terms of the hrst N orthogonal polynomials. This is 
very similar to (2.7). A crucial difference though between (2.7) and the more familiar expression 
(2.9) is that, in (2.7), the kernel does not depend on N, and the resulting matrix models are 
quite special. Matrix models of the form (2.7), coming from integral kernels of operators, have 
been considered before in for example [28, 29], and they have been recently studied from a more 
general point of view in [30], where the notion of “M-theoretic matrix model” was introduced. 

As we will see in the next sections, in the case of trace class operators obtained by quantiza¬ 
tion of mirror curves, the matrix model (2.7) is closely related to Chern-Simons matrix models 
[2] and to the matrix models appearing in the localization of Chern-Simons-matter theories 
[10, 31, 32] (see [12] for a review). In particular, the resulting matrix models are M-theoretic, 
in the sense of [30]; they have a well-defined ’t Hooft expansion in the limit (1.1), but they also 
have a well-defined M-theory limit, in which N is large but h is hxed. 

2.2 Operators from mirror curves and matrix models 

In [5], building on previous work on the quantization of mirror curves, it was postulated that, 
given a mirror curve to a toric CY threefold, one can quantize it to obtain a trace class operator 
p. As in [5, 14], we will focus on toric (almost) del Pezzo CY threefolds, which are defined as the 
total space of the anti-canonical bundle on a toric (almost) del Pezzo surface S, 

X = 0{-Ks) S. (2.10) 

The mirror curve depends on k complex moduli Za, a = 1, ■ ■ ■ ,k and is of genus one. This means 
that the complex moduli involve a “true” geometric modulus u and a set of “mass” parameters 
i = 1, • • • , r, where r depends on the geometry under consideration [33, 34]. The mirror curves 
can be put in the form 

IT(e^e^) = Os(x,y)+u = 0, (2.11) 

where Os{x,y) has the form 


k+2 

Os{x, y) = ^ exp (u["''x + v^^'^y + /i(Ci)) , (2.12) 

i=l 

and are suitable functions of the parameters The vectors ean be obtained from the 

fan defining the toric CY threefold. 
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In this paper we will be particularly interested in local where the function Os{x,y) is 
given by, 

C)p 2 (x, 2 /) = e*+ e^+ (2-13) 

and local F 2 , where it is given by, 

Of, (x, y) = e" + e^ + + ^e"". (2.14) 

To quantize the mirror curve (2.11), we promote x, y to self-adjoint Heisenberg operators x, y 
satisfying the commutation relation 

[x,y] = \h. (2.15) 

This promotes Os{x,y) to an operator, which will be denoted by O 5 (possible ordering ambi¬ 
guities are resolved by requiring the resulting operator to be self-adjoint). As conjectured in [5] 
and proved in [14] for many geometries, the inverse operator 

PS = O 5 ' (2.16) 

is positive-definite and of trace class. By the construction explained in the previous section, the 
fermionic traces of this operator, which we will denote by Zs{N,h), have an integral representa¬ 
tion in terms of the kernel of the operator ps- Therefore, in order to write down the matrix model 
(2.7), we need an explicit expression for this kernel. In general, obtaining such an expression is 
not easy. However, in [14], this problem was solved for three-term operators of the form 

0^,„ = e>^ + ey + e-“A (2.17) 


Geometrically, this operator corresponds to the anti-canonical bundle of the weigthed projective 
space F{1, n,m). When n = 1, these geometries can be constructed as partial blow-ups of the 
orbifold C^/Zm+ 2 ; where the orbifold action has the weights (1,1,m) [35]. In principle, more 
general geometries can be obtained by perturbing the operator Om,n appropriately. Note that 
the case m = n = 1 gives local P^, while m = 2 , n = 1 gives the ^ = 0 limit of local F 2 , i.e. a 
partial blowdown of the local F 2 geometry. Let us define 


Pm,n 


0 


-1 


(2.18) 


As shown in [14], the kernel of pm,n involves in a crucial way Faddeev’s quantum dilogarithm 
<kb(x) [36, 37] (in this paper, we use the notations of [14] for this function). We will also need 
the function 


^27Tax 

^a,c{x) 


(2.19) 


which behaves at large |x| as 


^a,c(a;) 


^—27TCX 


^27Tax 


X 00 , 

X — 00 . 


( 2 . 20 ) 


In [14] it was shown that, in terms of an appropriate variable p, related to x, y by a linear 
canonical transformation, one has 


Pm,n{pijP2) — (pi\Pm,n\P2) 


_ ^a,e(Pl) ^a,c(P2) _ 

2 bcosh(f(yi - P 2 ) -b 'f{a + c-nc)) 


( 2 . 21 ) 
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In this equation, the parameter b is related to h as 


,2 {m + n + l)h 

b =-, 

2tt 


while a, c are given by 


mb 


a = 


2{m + n + 1) ’ 


c = 


2 (m + n + 1)' 


( 2 . 22 ) 


(2.23) 


Since we have an explicit formula for the kernel of pm,n, we can write down an explicit expres¬ 
sion for the integral calculating the fermionic trace of this operator, which we will denote as 
Zm,n{N, h). This expression can be put in a very convenient form if we use Cauchy’s identity, as 
in [38, 39], 


n 


i<j 


2 sinh 


2 


2 sinh 




- = det- 


1 


2 cosh 

E n 


(2.24) 


cSSiv 


2 cosh 


2 


In this way one obtains, 

1 


^m,n(iV,h) = 


N\ 


where 


aN ^ 
2=1 


Cm.n. — 


^|2 ni<i4sinh(f(p^ -Pj))^ 

HiJ 2 cosh (f (pi - Pj) + \7rCm,n) ’ 

m — n -|- I 


(2.25) 


(2.26) 


2 (m + n + 1)' 

The above integral is real, since the kernel (2.21) is Hermitian. Although it is not manifest, it 
can be also checked that it is symmetric under the exchange m <-)• n. 

We are now interested in studying the matrix integral (2.25) in the ’t Hooft limit (1.1). 
Therefore, we should understand what happens to the integrand in (2.25) when h (or equivalently 
b) is large. To do this, we first change variables to 


and introduce the parameter 


27r 

Ui = -r-Pi, 


1 m + n + 1 1 




(2.27) 


(2.28) 


27r b^' 

Note that the strong coupling regime of h is the weak coupling regime of g. The crucial property 
to understand this regime is the self-duality of Faddeev’s quantum dilogarithm. 


4>b(x) = 4>i/b(x). 


(2.29) 


Then, we can write 


4/ + + <l>i/b((^ + 27ri(a + c)/b)/27rb 

"^a,c{P)\ e ^ I 9 -,rrr / <T). _L /ki /O-Trh-li ’ 


\27rgJ $i/b((tt-27ri(a-bc)/b)/27rb-i) 


30) 
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Figure 2: On the left, we plot the potential (2.34) for m = n = 1 (relevant for local P^), while on the 
right we plot it for to = 2 , n = 1 (relevant for local F 2 ). 


where u and p are related through (2.27). When b is large, 1/b is small and we can use the 
asymptotic expansion (see [40] for this and other properties of the quantum dilogarithm), 




(2.31) 


where is the Bernoulli polynomial. We dehne the potential of the matrix model as, 

Vm,n{u,g) = -glog|4'„,c(p)P, (2.32) 

where u and p are related as in (2.27). By using (2.31), we deduce that this potential has an 
asymptotic expansion at small g, of the form 


Vm,n{u,g) = 


t>0 


The leading contribution as g —)> 0 is given by the “classical” potential. 


27r 27r2 

By using the asymptotics of the dilogarithm, 

—x^/ 2 , X —)> 00, 


we find that 


Li2(-e") 




—e 


X —)■ —00, 


00 , 

-^U, u^-oo. 


(2.33) 


(2.34) 


(2.35) 


(2.36) 


Therefore, this is a linearly confining potential at inhnity, similar to the potentials appearing in 
matrix models for Chern-Simons-matter theories [39, 30]. The classical potentials for the cases 
m = n = 1 (relevant for local P^) and for m = 2 , n = 1 (relevant for local F 2 ) are shown in 
Fig. 2. We can now write the matrix integral as 


N 




m Un (27r)^ f \ 


ni<j4sinh(^^lV^) 

Hi 3 2 cosh (+ iTrCr, 


(2.37) 
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In the regime in which h is large (in particular, in the’t Hooft limit), we can use the asymptotic 
expansion (2.33) of the potential to study this matrix integral. The expression (2.37), computing 
the fermionic traces of the operators pm,n, is very similar to matrix models that have been studied 
before in the literature. The interaction between eigenvalues is identical to the one appearing in 
the generalized 0(2) models appearing in [41], and in some matrix models for Chern-Simons- 
matter theories studied in for example [38]. The parameter g corresponds to the string coupling 
constant, and the potential depends itself on g. Note however that, in the planar limit, only the 
classical part of the potential (2.34) contributes^. 

2.3 Perturbative expansion 

The matrix model (2.37) admits a standard’t Hooft expansion, of the form 

= \ogZ^^n{N,H) = ^h2-29j-M(A), A = ^. (2.38) 

9>0 

This can be easily seen by noting that, if we just keep the classical part of the potential, we have 
a generalized 0(2) matrix model, of the type considered in [41]. The corrections to the potential 
involve even powers of g, therefore they lead to corrections which preserve the form of (2.38). 

We would like to compute the genus g free energies Ideally, one would like to 

obtain them in closed form, as functions of the ’t Hooft parameter A, and this might feasible 
via a suitable generalization of the techniques of [41]. In this paper we are interested in testing 
whether the above free energies reproduce the genus expansion of the topological string, and we 
will perform a more pedestrian calculation, by doing perturbation theory in g to the very first 
orders. This means that we regard (2.37) as a Gaussian Hermitian matrix model, perturbed by 
single and double trace operators. The resulting perturbative expansion can then be converted, 
by standard means, into a weak coupling expansion around A = 0 of the very first 
This is very similar to the calculations done in [4] for the lens space matrix model. 

In order to work out the expansion of the matrix model, we first have to expand the “classical” 
potential Vm,h{u) around its minimum. Let us introduce the parameter 

9 = exp r — ) . (2.39) 

\m + n + 1/ 

Then, the minimum of the classical potential occurs at 

u* = '^ogXmiq), (2.40) 

where 

m _ -m 

Xm{q) = - —■ (2.41) 

q-q ^ 

The value of the potential at the minimum is given by 

FmlM = - ^logXm - ^ ImLi2 • (2.42) 

®The potentials appearing in some of the matrix models considered in [22, 23, 24, 25, 26] also involve (quantum) 
dilogarithms, but this does not necessarily indicate a deep connection between these matrix models and ours, due 
to the reasons listed in the introduction. 
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This can be written in a more compact form by using the Bloch-Wigner function 


D{z) = ImLi 2 (z) + arg(l — z) log \z\, 


(2.43) 


where arg denotes the branch of the argument between —tt and vr. One finds, 

(2.44) 

One can use the properties of the Bloch-Wigner function 

D{z) = -D{z), D{1- z) = -D{z), (2.45) 


to verify that (2.42) is symmetric under the exchange of m and n. For example, we have 

y®(ii* = 0) = |^, (2.46) 

where 

V = 2lm (hh (2.47) 

is related to the volume of the figure-eight knot. Similarly, we find 

^2,? = \ log(2)^ = (2.48) 

where C = 0.915966... is the Catalan number. As we will see, in order for the conjecture of [5] 
to be true, (2.42) has to be related in a precise way to a natural constant appearing in special 
geometry, namely the value of the large radius Kahler parameter at the conifold point. 

We can now expand the full potential appearing in (2.37) around u = «*, and write the 
interaction term as a “deformed” Vandermonde term. 


ni<j4sinh(^^^) 

Oij 2 cosh + iTrCm,n'j 


4:11 


-A\ui) (1 + 0 (A2(i,i))) , (2.49) 


where 

- Uj 

i<j 

is the usual squared Vandermonde, and 


(2.50) 


Im,n = COs(7rC'm,n) = sin 


vrn 


m + 1 + n 


(2.51) 


At leading order, we find a Hermitian Gaussian matrix model, and by including the corrections 
coming from the deformed Vandermonde and the potential (2.33), we can compute systematically 
the ’t Hooft expansion of Zm,n around A = 0. By using the standard formula for the partition 
function of the Gaussian matrix model. 


ZG{N,g,) = 


d 


N 




N\J (27r) 


N 


YliVi - Vjf = 


i<j 


N ‘^12 

9s 

(27r)^/2 


G{N + 1), 


(2.52) 
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where G{z) is Barnes’ function, as well as its asymptotic expansion at large N, we obtain the 
following results for the expansion of the genus g free energies appearing in (2.38). For 

the planar free energy, we find 




^flog__ 

2 V {m + n + l)I^^^A2 2) 


k=3 


(2.53) 


In this equation, Im,n is given in (2.51), A 2 is given by 


A 2 


sin 


(- 


r') sin (- 




27r- 


\^m+n+l J \^m+n+l J 


Sin 


(- 




\^m+n+l J 


(2.54) 


and the values of the coefficients /o,a: can be calculated explicitly as functions of m, n. The results 
for the very first k can be found in Appendix A. Similarly, one finds 


- 00 

^(---)(A) = -- log(A;i) + C'(-l) + £ /i,feA' 


j-(m,n)(A) ^ 


B 


k=l 


(2.55) 


29 


2g{2g - 2 ) 


^2-2g^^f^kX\ g>2. 


k=l 


The values of fg^k for the very first g, k, for general m, n, can be also found in Appendix A. We 
now list some results in the case of m = n = 1, relevant as we will see for local P^. We find, 


5671*’ 


^ flog f - ^A3 + ^A" + ^ A^ 

° 2 V \9V3 J 2 ) 47r2 9^/3 486 10935^/3 

10587r® 


A‘’ + O(A0, 

492075 ^ ^ 

Aa) = - A iog(A.)+c'(-i) 




(2.56) 


(M)ma_ 1 A-2 , 4776 ^ 3 1 8771 %2 


^2 (-^) 240^ 405%3^ 492075 


A^ + 0 (A^). 


For m = 2, n = 1, which is relevant for local F 2 , we obtain, 

r{ 2 ,l) 


A 2 

2 


tt^A 

4 


2(7, . 3 4 771® 5 73371® 


-^^’"(A) = w {log { ^ ) - h ) - ^A - -A® + :%^A^ - ^A® + V^A® + O(A^), 


71^ 

^2 


2 

A4 _ 

12" ' 288 960'' ' 172800 

T'Aa) = -4 log(A.) + C'(-l) + % - 

4"7a) = -Ta- + |4a-|^a^ + 0(a3). 


(2.57) 
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3. Testing the matrix model 


3.1 What should we expect? 

The conjecture of [5] gives a very precise prediction for the ’t Hooft expansion (1.2) of the 
matrix models arising from the trace class operators. To understand this prediction, we have 
to summarize some of the results of [5]. According to the conjecture of [5], the basic quantity 
determining the spectral properties of the operator ps is the modified grand potential fi). 

This grand potential depends on the “fugacity” p, which is related to the variable n entering in 
(2.5) as 

K = e^, (3.1) 

as well as on the parameters appearing in the operator Os, which we will collect in a vector 
The modified grand potential is determined by the enumerative geometry of the CalabmYau 
X. We first need a dictionary between the parameters p, and the parameters appearing in 
the enumerative geometry of X. First of all, we remember that the mirror curve (2.11) involves 
a modulus u. Typically, in mirror symmetry one uses the related modulus 


z = 


w 


(3.2) 


where, the value of r is determined by the geometry of X. There is a corresponding flat coordinate 
T, determined by 

—T = \ogz + 0{z). (3.3) 

Let us now consider the Kahler parameters of a general local del Pezzo, Ti, in an arbitrary basis. 
They can be always be written as a linear combination of the Kahler parameter T corresponding 
to z, and the Kahler parameters Tg., which are algebraic functions of the parameters appearing 
in the mirror curve. We have then. 


Tt^‘^T + a.,p. (3.4) 

The dictionary is now given as follows. On top of the standard mirror map (3.3), there is a 
quantum mirror map [7] of the form 

-T{h) = log( 2 :) + '^ai{h)z^. (3.5) 

£>1 

The “effective” p parameter is defined by an expansion similar to this one, 

PeS = p-- (3.6) 

The Kahler parameters Tj are then related to the parameters p, by the following equation, 


Ti = CiPes + (XijT^.. (3.7) 

We can now write down the general expression for the modified grand potential. It is of the form, 

K) = J^'^\peS,ii, K) + JM2{PeS,$,, K) + JwS(l^eff, 4, ^)• (3.8) 
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(3.9) 


Here, is the perturbative part, which is a cubic polynomial in 

h) = ^ ^ + Bi{^)h}j n + A{^, h). 

The constants appearing in this expression can be obtained from a semiclassical analysis of the 
operator ps (see [5] for details). The function A{^,h) is not known in closed form, but in some 
examples there are educated guesses for it. The “membrane” part of the potential has the form, 

'^M2(/^eff) 4) ^) heff'^b(/^eff) ^) ^) T </c(A*eff) ^) ^)) (3.10) 


where 


Jbil-^eS} fi) — 2 ^ E E ^ d) 


sinf^(2jL + l)sinf (2j^ + l) 

6 1 


Jl^r w,d 


2w sin^ ^ 


7/ c ^ Ard sin^(2ji + l)sin^(2jK + l) _^d.T . 

Mfies, = - 2tcsin3^ - 


i,j JlJr u^.d 




27]- ^ " Qfi 

JlJr ui,d 


sin^(2ji + l)sin^(2j^j + l) 


hw I 

_^ 

2hw^ sin3 ^ 


jyd -U.dT 
JL,JR^ 


In this equation, are the rehned BPS invariants of the CY Y, T = (Ti,T 2 ,---) is the 

vector of Kahler parameters, c is the vector of constants Cj appearing in (3.4), and d is the 
vector of degrees. Finally, the worldsheet part of the modified grand potential is 

-''ws(/refr,C,h) = r2sin?^'j (-l)‘^'^e“^"‘^''^. (3.12) 

g>0 d,v ^ ^ d 


In this equation, are the Gopakumar-Vafa invariants of Y, while B is a B-field, which is 
related to the anticanonical class of S', and necessary for the cancellation of poles in Js{p,$,h) 
[42]. 

According to the conjecture of [5], the spectral determinant of the operator ps is given by 

Es{p, h) = Y^ e'^s(/^+27rm4,;i)_ 

nGZ 


As first shown in [43], in the context of ABJM theory, this representation leads to a very con¬ 
venient formula for the fermionic trace Zs{N,h) as an integral transform of the modified grand 
potential, 


Zs{N,h) 


1 

27ri 




(3.14) 


The contour C appearing in this integral is shown in Fig. 3, and in view of the cubic behavior of 
JsifJ', K), it leads to a convergent integral (this is the contour used to define the Airy function). 

As already pointed out in [5], the above conjectural results lead to a very precise formula 
for the’t Hooft expansion of Zs{N,h). To obtain this expansion, we follow a procedure similar 
to what was done in [44]. We first note that the function Js{p,$,h) has itself a ’t Hooft limit in 
which 


h 

h 


/U —)• oo 


h —)• oo 


C fixed. 


(3.15) 




Figure 3: The contour C in the complex /x plane, which can be used to calculate the fermionic trace 
Z{N,h) from the modified grand potential. 


It is clear that this double scaling is needed if we want the integral in the r.h.s. of (3.14) to be 
non-trivial as /i —)■ oo. What is the effect of this limit on ^)? First of all, we have that 


h 


C, 


(3.16) 


and all the exponential corrections to (3.6) vanish. Similarly, the membrane corrections in (3.10) 
vanish as well, since they are exponentially small as ^ oo. The only surviving terms come 
from the perturbative and the worldsheet parts of the modified grand potential. As we will see 
in examples, the function A(^, h) has an expansion as li —)■ oo of the form 


= (3-17) 

g>o 

where Ai(^) includes as well a logarithmic dependence on h. We conclude that, in the ’t Hooft 
limit (3.15), the modified grand potential has the expansion, 

OO 

= (3.18) 

3=0 


where 


JgiCX) 


bvr 


^3 + BiC + Ao(^) + 


D{^) 

h 


+ 


^ Z 7 inst 

IbTT^^O 



BoC + A,{^) + Fr^(^t,^y 

A,(0 + (47r2)2s-2i.inst 


Here, we have introduced the variable 

t = 2TTr(, 


(3.19) 


(3.20) 


and is the standard genus g topological string free energy, expressed in terms of 

the conventional Kahler parameters t, (these correspond to the parameters T, T^,, up to a 
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rescaling by ‘1'k jK). We now note that, in order to perform the ’t Hooft limit, we must also make 
a choice for the scaling of the parameters T^.. If these scale with fi when /i —)• oo, the ’t Hooft 
limit of Jg{C:0 reproduces the standard topological string free energy at genus g. However, we 
might also want to keep the T^. fixed. In that case, in order to obtain the large h expansion, we 
have to re-expand the functions Jg{C,0^ one finds 

OO 

= (3.21) 

3=0 


where, for the very first orders, 


Jlic.i) 




+ i^rr(t,o). 


(C, = BoC + Ai(^) + F'r^ (t, 0) + 


a^Fr^(t,o) 

87r2 


(3.22) 


Similar expressions can be of course obtained for the higher order functions Jg{Cj^): involving 
derivatives of the with g' < g. Here we have assumed that derivatives of the free energies 
w.r.t. an odd number of parameters and evaluated at = 0, vanish, in order to have an 
expansion involving only even powers of h. This is the case in the example of local F2, which 
will be analyzed in detail later on, but might not be a universal fact. The formulae above can 
be easily modified if this is not the case, and the conjecture of [5] would then imply that the 
’t Hooft expansion of the fermionic traces has odd powers of h. This does not happen for the 
operators pm,n that we are analyzing in this paper, but might happen in other situations. 

In any case, the formulae of [5] lead to precise expressions for the ’t Hooft limit of the 
modified grand potential. One interesting aspect of this limit is that it only keeps the worldsheet 
instanton corrections, which are precisely the most difficult ingredients from the point of view of 
spectral theory. In order to obtain the ’t Hooft expansion of the fermionic trace. 


log Zs{N,h) = (3.23) 

3>0 


we still have to calculate the integral in (3.14). It can be evaluated in a systematic asymptotic 
expansion at large h by using the saddle-point approximation. At leading order, we find that the 
’t Hooft parameter is given by 

dJl^ 

dc 

This determines A as a function of C, and conversely, () as a function of A. The genus zero free 
energy is then given, as usual, by a Legendre transform, 




(3.25) 


We obtain in particular 

The next-to-leading order correction T'f^(A,^) is given by the one-loop approximation to the 
integral (3.14), 

.^f(A,€) = jf(C(A),0 - ^log (2vr^) . (3.27) 
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Higher corrections can be computed in a straightforward way, and we note that, when the 
parameters T^. are fixed in the’t Hooft limit, we have to use the modified formulae (3.22). 

The calculation of higher corrections is better understood if we realize that the integral (3.14) 
has precisely the form appropriate to implement a symplectic transformation of the topological 
string free energies, as explained in [45]^. The equations (3.24) and (3.26) indicate that the 
integral (3.14) is implementing an S transformation, in which the variable C in becomes 

the dual variable From the point of view of the topological string, the 

are free energies in the so-called large radius frame (where we have an interpretation in terms of 
a worldsheet instanton expansion). The S transformation, in these local del Pezzo geometries, 
takes the large radius frame to the conifold frame. We conclude that, according to the conjecture 
of [5], the functions T'^(A,^) appearing in the ’t Hooft expansion of the fermionic traces (i.e. in 
the ’t Hooft expansion of the matrix model (2.7)) are the genus g free energies of the topological 
string in the conifold frame (up to normalizations and integration constants, as we will see). 

The appearance of the conifold frame is completely natural in view of the matrix model 
representation. Topological string theory in the conifold frame is known to have a universal 
structure: when the free energies are expanded around the conifold point, the leading singularities 
are of the form tc ^, where tc is the flat coordinate around the conifold, and the coefficients 
of these singularities are determined by the c = 1 string free energy [46]. There is in addition 
a “gap” condition, [47], which says that the corrections to this leading singularity involve only 
non-negative powers of tc- It has been known for some time that this structure is precisely what 
is found in the weak ’t Hooft coupling expansion of a perturbed Gaussian Hermitian matrix 
model. Therefore, the appearance of the conifold frame in the above derivation is already what 
one should expect from the existence of an explicit matrix model description. 

In summary, if the conjecture of [5] is true, the matrix model (2.7) has a’t Hooft expansion, 
given by the genus g topological string free energies in the conifold frame. We will now test this 
prediction in detail in the case of the matrix integral (2.37). For n = m = 1, this is just local 
while for n = 1, m = 2 this is local F 2 for a specific value of the “mass” parameter, namely 
^ = 0. In other words, according to the conjecture of [5], we should have 

T-(i’i)(A) = J-f (A), T-f')(A) = (A,^ = 0), (3.28) 

for all g > 0. The l.h.s. of these conjectured equalities has been computed, in an expansion 
around A = 0, in (2.56) and (2.57). In the next two sections, we will use topological string theory 
to compute the r.h.s. of (3.28) and check to the available order that these two functions are 
indeed equal. 

3.2 Local 

The modified grand potential of local P^ has been determined in [5] in great detail. In this case, 
there are no mass parameters In order to write down the results, let us recall some elementary 
facts about the special geometry of local P^, relevant for our calculation. The Picard-Fuchs 
equation determining the periods is 

(03-3z(30 + 2)(30 + l)0)n = O, (3.29) 

where 

e = z-^, (3.30) 

dz 

^Note however that the integral transform appearing in [45] is purely formal, and relates two asymptotic power 
series, while the integral (3.14) has a precise non-perturbative meaning. 
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large radius 
z = 0 



conifold 



Figure 4: The moduli space of local can be parametrized by a single complex variable z, and it has 
three special points: the large radius point at z = 0, the conifold point at z = 1/27, and the orbifold point 
at z = oo. 


and z parametrizes the moduli space, which has three special points: z = 0 is the large radius 
point, z = 1/27 is the conifold point, and z = oo is the orbifold point (see Fig. 4). A basis of 
solutions around the large radius point z = 0 is given by 


a;i(z) = log(z) + 6z + 45z^ + 560z^ + C)(z^), 

(6. + 45.= + 560.= + 0(.‘‘)) + 3 . + 

DO 4 d 

These periods determine the large radius, genus zero free energy Fo(t) as 

-t = iOl{z), 

dFo 


(3.31) 


dt 


= (^2{z), 


which gives 


m) = ^ - 3e-* - ^e-2* - + ©(e"^^) 


(3.32) 


(3.33) 


Note that the signs are not the standard ones. This is due to a non-trivial i3-field which has to 
be turned on in (3.12) in order to obtain a consistent modified grand potential, as first noted in 
[42]. 

We can now write down the ’t Hooft limit of the modified grand potential. The constants 
appearing in the perturbative piece (3.9) are [5] 


2 ’ 


= 0 , 


Bo = 

u 2’ 


Bi = - 

Ibvr 


On the other hand, the function A{h) is given by 

3Ac{h/Tr) — Ac{3h/Tr) 


m = 


(3.34) 


(3.35) 


where [48] (see also [49]) 

Ac{k) = 


2C(3) 

Ti'^k 


16y vr^ Jq e^* — 1 


log(l — e ^)dx 


(3.36) 
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Figure 5: The’t Hooft parameter X = N/h a.s a. function of the modulus 0 < z < XjTl. The large radius 
point z = 0 corresponds to strong ’t Hooft coupling, while z = 1/27, the conifold point, is the weakly 
coupled theory. 


This function has a large k expansion of the form, 

= -^C(3) + - log(2) + 2C'(-1) + g log (-) + ^ - (3.37) 

^ S>2 ^ ^ 


where 


We conclude that 


c„ = 


B2gB2g-2 


{4g){2g-2){2g-2)V 


(3.38) 


At) = 


3C(3) ^ 


^1 = log(^) + C'(-l) + ^ log(27r) + ^ log(3), 
Ag = (47r2)29-2 (3 _ 32-23) {-l)9-^cg, g>2, 


(3.39) 


The function Jq^ {C) is expressed more conveniently in terms of the variable t appearing in (3.20), 
which is now given by, 

t = OttC, (3.40) 

since r = 3 in this geometry. We find. 








Fo(t) + 3C(3) — 


(3.41) 


The’t Hooft parameter is now given by (3.24), which reads in this case, 

StT^ _ dFo TT^ 

~ 1h~ Y' 


(3.42) 


The r.h.s. of this equation is nothing but the vanishing period at the conifold point. Therefore, 
the ’t Hooft parameter varies between 0 and 00 as z varies between 1/27 and 0, as shown in 
Fig. 5. The region around the large radius point z = 0 in CY moduli space corresponds to 
strong’t Hooft coupling, while the region around z = 1/27, the conifold point, corresponds to 


- 18 - 








the weakly coupled theory. Since we are interested in expanding the free energies around A = 0, 
we have to analyze the theory around the conifold point. To do this, we define the variable 


y = l- 27z. 


(3.43) 


We now solve the Picard-Fuchs equation near the conifold point, i.e. near y = 0. There are 
again two independent periods. One of them, which we will denote by tc{y), is a flat coordinate 
near the conifold point. It is given by the power series expansion 


lly2 I09y3 9389/ 88351/ 823187/ , 7 . 

^c{y) yT .|o T o/to oen a a oncn /i c oa ooc A c V^/’ 


18 243 26244 

and is related to the’t Hooft parameter as 

/3 


295245 3188646 


A = 


127r2 


tc{y)- 


(3.44) 


(3.45) 


Therefore, as announced above, the ’t Hooft parameter defined by the fermionic traces is a flat 
coordinate at the conifold point, proportional to tc{y)- 

The period uJi{z) defines the genus zero free energy T'g^(A). Indeed, from (3.26), (3.40) and 
(3.32), we find 

// . (3.46) 

dX dvr 

This function is, up to normalizations and integration constants, the genus zero free energy of 
local in the conifold frame, and it has been computed in for example [50]. In order to expand 
it around A = 0, we need the expansion of the period Ldi{z) around the conifold point. This is a 
standard exercise in special geometry and one finds. 


\/3 / 

ioi{z) = + — f tc(y) log 


+ s (y) > 


where 


s(y) ^2 1458 

and the constant c is given by 


4 5 


31og(3) + 1 

7y2 877/ 176015/ 9065753/ 17960917/ 


(3.47) 


+ 


314928 


+ 


17714700 


+ 


38263752 


+ O (/) , (3.48) 


9 

2.90759 


V3. 


y2,2 


c = -- 4 T 3 ( 1,1, -, 2, 2, 2; 1 ) + 3 log(3) = Re { ( -1 


27r 


1 2 ^ 

3’ 3’ 

0 , 0,0 


(3.49) 


The relationship (3.46) can be integrated to obtain the genus zero free energy, up to an integration 
constant. This constant can be determined as follows. The value of J-q{X) at A = 0 is given by 


(0) = 4 (C(0)) = 


1 




TT^C 


^o(c)- —+ 3C(3)), 


where we used that 


t{z = -]=c. 


(3.50) 


(3.51) 
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The value of Jq^(C( 0)) can be obtained numerically with high precision and we find that J-'q^(O) 
vanishes^ i.e. we find that 


i^o(c)-^ = -3C(3). 
When all this is taken into account, we find the expansion, 


(A) 


c 2 /log(A) 3 f 2 tt W tt^A^ tt^A^ Sbvr^A® 
^ 4 ^ V^y y ~ 9v/3 ^ 486 ^ 10935\/3 

IOSStt^A® 33927r^°A^ 2087447ri2;^® . 9^ 

492075 2066715\/3 ~ 1171827405 ^ 


(3.52) 


(3.53) 


This expansion should exactly agree with the expansion of calculated in (2.56), and it 

does, provided the number c given in (3.49) satisfies 


c = 


9V 


(3.54) 


Amazingly, this is a true identity [51, 52]! (see also [53]). This encapsulates the power of the 
conjecture of [5]; the constant c comes from the world of topological string theory, and it gives 
the value of the Kahler parameter t at the conifold point. The constant V comes from the world 
of trace class operators arising from mirror curves, which lead to quantum dilogarithms by the 
results of [14]. For the conjecture of [5] to work, the two numbers coming from these two different 
worlds have to agree. And they do. 

The agreement between T'q^(A) and Tq^A)(_\) (at the order we have been able to check) is a 
non-trivial analytic test of the conjecture put forward in [5]. Specifically, this calculation tests 
one of the most surprising claims of [5], namely that the non-perturbative corrections appearing 
in the spectral theory of the operator pp 2 are encoded in the standard topological string. 

Let us now consider the genus one free energy. By using well-known results for the B-model 
of local P^, one finds 


«^r(C) = ^ log log (1 - 27z)) + Ai. 


Using (3.27), as well as 


,2 Tp2 


d^J 


-dtr 




(3.55) 


(3.56) 


one obtains the following expansion for IFf (A), 

p2 1 , bvr^A tt^A^ 407r®A^ 2837r®A^ 13767r^°A^ 

1 ( )- ;l2 ”^18^/3 486 2187^/3”^ 32805 177147^/3 

72272tt^‘^X^ 79367ri^A^ . g. 

^ 55801305 ^ 14348907v/3 ^ ^ 


(3.57) 


which again is in precise agreement with what was found in (2.56). 

There is a further test of the first conjectural equality in (3.28) which can be made at higher 
genus, by using the results of [50]. Let us denote by Fg(t) the higher genus free energy of local 
P^, as computed in [50], with an extra sign (—1)^“^. Then, one has 

= g>2. (3.58) 
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The genus g free energy appearing here, Fg{t), includes the so-called constant map contribution, 

Fgit) = 3{-iy-^Cg + O (e-*) . (3.59) 

In our formalism, this contribution comes from the first term in the r.h.s. of (3.39), while the 
second term leads to an additional constant in (3.58). As mentioned above, the higher genus 
functions Fg (A) should be given by the symplectic transformation of the (3.58) to the conifold 
frame. This was done for Fg{t) in [50]. The resulting quantities, when expanded around the 
conifold point, display the singular term in appearing in (2.55), plus a constant, and a 

series starting in A, i.e. 

= 2a{2g- + 0((,). (3.60) 

Here, tc is a local flat coordinate around the conifold, proportional to A. The constant term 
appearing here cancels exactly the second term in the r.h.s. of (3.58). This is precisely what is 
required by (2.55). In other words, the function A{H), which was conjectured in (3.35) in order 
to reproduce the spectral properties of /9p2, is precisely what is needed to guarantee the matrix 
model behavior of (A) near A = 0. It is also easy to verify that, after taking into account 

the appropriate normalizations, the expansion of T'2^’^^(A) agrees with the genus two topological 
string free energy at the conifold point. 

As a final remark, note that the planar free energy T'q^(A) can be also computed at strong 
’t Hooft coupling, and an elementary calculation gives 

Ff{\)^-^\^/\ A»l. (3.61) 

This is the typical behavior for a theory of M2 branes [54], and it agrees with the M-theory limit 
of the free energy computed in [5]. The behavior of the matrix model for local is therefore 
very similar to what is found in the ABJM matrix model and its generalizations [55, 56, 39]. 

3.3 Local F 2 

The case of local F2 is an interesting one, since we have an extra parameter ^ in the operator 
(2.14). Correspondingly, there is an extra Kahler parameter in the geometry, which is usually 
denoted by Tb, and related to by 


i = 2 cosh 


Tb 


(3.62) 


The perturbative part of the modified grand potential has been determined in [57] 




SttH 


F + 


I ( TT 


h 


-log 

3 27r ^ 


2 


h 


g + (3.63) 


although the function A{h,^) is not known in full generality. The genus g free energy of this 
geometry, t^); which appears in the expressions (3.19), is now a function of two variables. 

The variable t, defined in (3.20), is related to C by 

t = 47r(^, (3.64) 
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since r = 2 in this geometry. The variable is is related to Tb by a rescaling of 27r/h, and it 
defines a parameter m by the analogue of (3.62), 

m = 2cosh^^^. (3.65) 

The genus zero free energy FQ{t,tB) can be determined, up to integration constants, by the two 
equations [58] 

dt _ 2 f 16z 

dz TTZy^l - 4(2 + m)z \4:{2 + m)z-l 

d^Fp ^ _ 2 f 4{m + 2 )z-l ^ 

dzdt zy/l - 4(m - 2)z \4:{m - 2)z - 1) ' 

where z is the global coordinate corresponding to t. 

In order to test the second conjectural equality in (3.28), we have to consider the theory for 
^ = 0. The corresponding value of Tb is 


(3.66) 


Tb = ivri. 


(3.67) 


Therefore, the’t Hooft limit of the modified grand potential involves the functions written down 
in (3.22). Since we have to expand ^b) around tB = 0, we have to consider the topological 

string theory on local F2 for the value of the parameter m = 2. For this value, the theory is 
identical to the diagonal, local x geometry, and (3.66) simplifies to 


dz 

d^Fp 

dzdt 


- K (16z), 

TTZ 

—K (1 — I62;). 


(3.68) 


This can be integrated explicitly, and one finds the two periods, 


u:i{z) = log{z) + Az 4 F 3 ( 1,1, 2,2,2; I62; ) , 


i02{z) = I 


TT 


i i 1 
2 ’ 2 ’ 

0 , 0,0 


3 ’ 


(3.69) 


which solve the Picard-Fuchs equation 

{ 9 ^ - Aze{2e +1)2) n = 0. 


(3.70) 


The genus zero free energy depends now only on a single parameter t, and we will denote it 
Fo{t) = Fo{t,0). It is determined again by (3.32), and one finds the expansion, 

To(t) = |-4e-‘-2e-2*---- (3.71) 

Collecting the above results, one finds that the function ^^(C)C = 0) is given by 

^^(C,e = 0) = (^o(t) - ^ + IQTT^AoiC = 0)) . (3.72) 
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Figure 6: The’t Hooft parameter X = N/h a.s a. function of the modulus 0 < z < 1/16. The large radius 
point z = 0 corresponds to strong ’t Hooft coupling, while z = 1/16, the conifold point, is the weakly 
coupled theory. 


The equation (3.24) determining the’t Hooft parameter reads in this case. 




dFoit) 

dt 



(3.73) 


As in the case of local 
located at 


the r.h.s. of this equation is a vanishing period at the conifold point. 



(3.74) 


The ’t Hooft parameter varies between 0 and oo as z varies between 1/16 and 0, as shown in 
Fig. 6. As in the case of local we are interested in expanding the quantities at weak’t Hooft 
coupling, near the conifold point. We introduce the local variable 


y = 1 — 16z. 

There is a flat coordinate near the conifold point y = 0, with the expansion, 
^ ^ 5y\ 89y^ , 381y^ , 25609y5 ^ ^ 

ic(2/)-y+^ + ^ + ^^ + ^j^ + 0(y ). 


It is related to the’t Hooft parameter by 


A = 


tcjy) 


The period cui(z) has the following expression near the conifold point, 

8C tc{y) . f y\ 1 

uji{z) = -+-log + -s(y)> 

TT vr Mb/ TT 


where 


35r 2141y4 465061^ 


s(y) = —y -1-^—h 

" 16 288 12288 

-F2 


+ 


2457600 


+ o {y^). 


The genus zero free energy (A) is determined by the equation 


5J-^^(A) _ cui(z) 


dX 


dTT 


(3.75) 

(3.76) 

(3.77) 

(3.78) 

(3.79) 

(3.80) 
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This can be integrated to give, up to an integration constant 


4^(A,e = 0) = -^A + A 


2 C 


TT^ 


+ 


7337r®A6 

172800 


^ log(A) 

477riOA7 

16128 


TT 


- 7 + log n 


+ O (A®) . 


tt^A^ Stt^A^ 
12 288 


77r®A5 

960 


(3.81) 


(<2 1 

This expansion agrees with the result for Jq ’ '^(A), as obtained from the matrix model in (2.57). 
The coefficient in the linear term in A involves the Catalan number, due to (3.78), and this is 
precisely what is required to agree with the value of the potential ^/^^(tt*) at its minimum. Full 
agreement between both expressions requires the integration constant for T'q^(A) to vanish. We 
have verified numerically that this will be the case if 


AoiC = 0) 


5C(3) _ 

IOtt^ 


(3.82) 


It would be interesting to determine the function A(h,^) and check whether or not it satisfies 
(3.82)5. 

Note that ^-^^(A,.^ = 0) is not the genus zero free energy of local F 2 in the conifold frame 
and with m = 0. Rather, it corresponds to the value m = 2. The reason for this “mismatch” 
is that, in the conjecture of [5], the Kahler parameters corresponding to the mass parameters 
appear divided by h in the worldsheet instanton expansion. Therefore, when /i —)■ 00 , they are 
frozen to the values which correspond to txi = 0. 

Let us now work out the next-to-leading order in 1/h^, which is given by the second line of 
(3.22). We first note that 


5Fo 

= 0, 

d‘^Fo 


IdFo 

dts 

is=0 

dtl 

tB=0 

2 dm 


(3.83) 


This last derivative can be obtained, as a function on moduli space, by taking a further derivative 
w.r.t. t and using (3.66); 


dmdy 


m=2 


(1.1;2; ^) + ^A(l - vmy) - 1) 


(3.84) 


where K{y), E{y) are the elliptic integrals of the first and the second kind, respectively, and we 
have expressed it already in terms of the conifold variable. From the general formula in (3.22), 
one finds 

: A , rn /j. n\ ^ *9Fo | 


•^1 (C) = ^+^i + i"i(L0)- 


16 dm 


m=2 


where 


Fi{t,0) =--log 


(:^4^)-Llog(642(l-16z)) 


Using again (3.27), but this time for the functions Jq^(C), we find. 


AP(A,« = 0) = --l„g(A) + 2_- 288 


tt^A Stt^A^ tt^A^ 537r®A^ 
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(3.85) 

(3.86) 

+ O (A®) , (3.87) 


® While this paper was being typed, Y. Hatsuda found an ansatz for the function A{h, ^ = 0) [59] which indeed 
leads to (3.82). 
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up to an additive constant independent of A (but depending on h through Ai ). This again agrees 
with the expansion in (2.57). Note however that the function J-']^^(A,^ = 0) is not the genus one 
free energy of local F 2 at the conifold point and for m = 0 (or even m = 2), since it involves 
additional terms. This reflects the fact that the relation between the parameters appearing in the 
operator and the “mass” parameters appearing in the geometry depends on h. In any case, the 
result we have obtained is in perfect agreement with the conjecture of [5], and gives a non-trivial 
test of the way in which mass parameters are incorporated in that framework. 

Although (3.87) is not the standard genus one topological string free energy of local F 2 , 
this does not contradict our claim that one can obtain these free energies from the mirror curve 
operators and their fermionic traces. If we want to recover the standard dependence on the mass 
parameters, we have to take a sort of Veneziano limit, and scale the parameters appearing in 
the operators in an appropriate way, as h grows large. For example, if we consider the fermionic 
traces Z{N,^,h) of the operator in the ’t Hooft limit (1.1), and we choose the following 
scaling with h for the parameter appearing in the operator, 

^ ^ (^) ’ 

the conjecture of [5] says that the resulting 1/N expansion would be governed by the standard 
topological string free energy on local F 2 , with a generic Kahler parameter tB^- 

4. Conclusions and open problems 

In this paper we have proposed a new type of matrix models whose 1/N expansion gives the 
all-genus free energy of topological strings on toric CY threefolds. Our proposal is based on the 
conjecture of [5], and it is conceptually clear and elegant, as it can be seen in Fig. 1: given a 
mirror curve, its quantization leads to a trace class operator, as shown in [14]. The fermionic 
traces of these operators admit an integral representation, and this leads to the matrix models of 
this paper. This construction provides a non-perturbative completion of the standard topological 
string, as it was made clear already in [5]. However, the representation in terms of matrix models 
spelled out in detail in this paper is particularly appealing, since it involves the standard’t Hooft 
limit which underlies the string/gauge theory correspondence. 

In practice, in order to study these matrix models, one needs an explicit expression for the 
kernel of the corresponding trace class operator. This was achieved in [14] in some cases, and 
the resulting matrix models are labelled by two numbers m, n. They describe topological string 
theory on the anticanonical bundle of the weighted projective space F(l, m,n), which can be 
obtained as limits of well-known toric geometries. In particular, the case m = n = 1 gives an 
explicit matrix model for local F^. 

The fact that the’t Hooft expansion of these matrix models agrees with the topological string 
free energy is an analytic, non-trivial test of the conjecture of [5], which in addition probes the 
non-perturbative sector of the spectral problem. In this paper we have performed two detailed 
comparisons, involving local and local F 2 , and found a complete agreement between the matrix 
model free energies and the topological string free energies. Interestingly, our matrix models give 
the free energies in the conifold frame. It has been suspected for a long time that this is the 
natural frame for a matrix model representation of the topological string, due to the conifold 
behavior of the free energies and the gap condition [47, 50]. Our proposal realizes this idea in 

®This expectation has been checked in detail in [60], after the first version of this paper was posted. 
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a very concrete way. Indeed, it implies the gap condition for the CY threefolds that we have 
studied. The BKMP conjecture [61, 62] should be also a consequence of our proposal, since we 
give an explicit matrix model realization of the topological string free energies, and very likely 
our matrix models satisfy the topological recursion of [63] (or a variant thereof). 

The matrix models obtained in this paper are very similar to the ones appearing in the 
localization of Chern-Simons-matter theories. We have now a complete parallelism between the 
theory of Chern-Simons-matter matrix models and the theory of topological strings on toric 
CY threefolds: in both of them, the perturbative sector is computed by the ’t Hooft expansion 
of a matrix integral, but there are non-perturbative corrections to the ’t Hooft expansion. In 
the case of the matrix models for topological strings considered here, these non-perturbative 
corrections are explicitly known, since they can be obtained from the modified grand potential 
of the theory. They correspond to the “membrane” part of the modified grand potential (3.10), 
which is determined by the NS limit of refined topological string, and they are exponentially 
small, of the form, 

exp (;i5aTo^(A)) , (4.1) 

just like the non-perturbative corrections to the ABJM matrix model obtained in [64, 39]. Note 
that, since /i ~ this is indeed non-perturbative in the topological string coupling constant. 

Clearly, there are many avenues for future research. Our analysis of the matrix models (2.37) 
has been very elementary, and based on a perturbative expansion. It would be very interesting to 
solve these models exactly in the planar limit and beyond, by using for example the techniques of 
[41]. We should also obtain explicit representations for the kernels of other operators appearing in 
the quantization of mirror curves, in order to write down explicit matrix integrals in more general 
cases. Another interesting problem would be to use the topological string free energies and their 
trans-series extensions, as constructed in for example [65, 66] to reconstruct the fermionic traces 
and their matrix model representation via resurgent analysis. In this paper, as in [5], we have 
focused on mirror curves of genus one. It would be important to understand in more detail how 
to generalize this construction to higher genus mirror curves. Finally, it would be interesting 
to know if the matrix integrals we are writing have some gauge theory interpretation, or some 
underlying M2 brane interpretation, as speculated in [5]. 
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A. Results for the perturbative expansion 

In this Appendix we list some results for the coefficients fg^k appearing in the weak coupling 
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expansion of the free energies in (2.53) and (2.55), for general m, n. One finds, for genus zero. 
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Finally, for genus two, we obtain: 
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